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ABSTRACT 

We study the hardness-intensity correlation (HIC) in gamma-ray bursts (GRBs). In particular, we 
analyze the decay phase of pulse structures in their light curves. The study comprises a sample of 82 long 
pulses selected from 66 long bursts observed by the Burst And Transient Source Experiment (BATSE) 
on the Compton Gamma-Ray Observatory. We find that at least 57% of these pulses have HICs that 
can be well described by a power law. A number of the other cases can still be explained with the power 
law model if various limitations of the observations are taken into account. 
Q-i The distribution of the power law indices 7, obtained by modeling the HIC of pulses from different 

bursts, is broad with a mean of 1.9 and a standard deviation of 0.7. We also compare indices among 
pulses from the same bursts and find that their distribution is significantly narrower. The probability p 
of a random coincidence is shown to be very small (< 2 x 1CU 5 ). In most cases, the indices are equal to 
within the uncertainties. These results demand a physical model to be able to reproduce multiple pulses 
with similar characteristics for an individual burst, but with a large diversity for pulses from an ensemble 
of bursts. This is particularly relevant when comparing the external versus the internal models. 

In our analysis, we also use a new method for studying the hardness-intensity correlation, in which 
the intensity is represented by the peak value of the EFe spectrum, where E is the energy and Fe is 
the energy flux spectrum. We compare it to the traditional method in which the intensity over a finite 
energy range is used instead, which may be an incorrect measure of the bolometric intensity. This new 
method gives stronger correlations and is useful in the study of various aspects of the HIC. In particular, 
it produces a better agreement between indices of different pulses within the same burst. Also, we 
find that some pulses exhibit a track jump in their HICs, in which the correlation jumps between two 
power laws with the same index. We discuss the possibility that the track jump is caused by strongly 
overlapping pulses. Based on our findings, the constancy of the index is proposed to be used as a tool 
for pulse identification in overlapping pulses and examples of its application are given. 

Subject headings: gamma rays: bursts - gamma rays: observations - methods: data analysis 
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method, it is essential to find common properties among 
pulses. 

The hardness-fluence correlation was discussed first by 
Liang & Kargatis (1996) who described it as being an 
exponential decay of the spectral hardness as a function 
of the photon fluence. The exponential decay constant 
appeared to be invariant between pulses in some bursts, 
which led the authors to suggest that the pulses are cre- 
ated by a regenerative source rather than in a single catas- 
trophic event. However, Crider et al. (1998a) dismissed the 
apparent invariance as coincidental, and consistent with 
drawing values out of a narrow statistical distribution, 
combined with rather large uncertainties in the determi- 
nation of the exponential decay constant. 

The hardness-intensity correlation was discussed first by 
Golenetskii et al. (1983, hereafter G83). In the present 
work, we study the HIC for the decay phase of a sample of 
GRB pulses observed by the Burst and Transient Source 
Experiment (BATSE) on the Compton Gamma-Ray Ob- 
servatory (CGRO). We present an extensive comparison of 
the HIC behavior among pulses from the same burst and 
between bursts, partly motivated by the behavior of the 
Hardness-Fluence Correlation. Some of the results have 
been presented in preliminary form in Borgonovo & Ryde 
(2000| and Ryde, Borgonovo, & Svensson (2000). 

First, in §0 we discuss previous work and results, in 
which the HIC was often found to be a power law relation. 
In §y, we present the data we used in the analysis and dis- 
cuss the observations (§3.1), the sample selection (§3.2), 
the spectral modeling (§3.3) and the analysis method used 
(§3.4). In particular, in §3.4.2, we introduce a new method 
to analyze the HIC. We present our results in §4. The use- 
fullness of our new analysis method is shown in §4.1, and in 
§4.2 we study individual pulses in GRBs and present the 
general distribution of their power law indices in §4.2.1. 
The cases that were not included in the analysis of the 
distribution, are examined in §4.2.2. In §4.3 we turn to 
the study of multi-pulse GRBs and investigate cases with 
several well separated pulses (§4.3.1) and discuss charac- 
teristic cases exhibiting track jumps in their HICs. In par- 
ticular, §4.3.2 presents cases with track jumps occurring in 
apparently single pulses. In §4.4 we study the power law 
HICs of pulses within the same burst and compare them 
to the general distribution from §4.2.1 and find that they 
are more alike than what is expected from the general dis- 
tribution. We discuss the final results of our analysis in §pl 
Our new analysis method is discussed in §5.1 and the track 
jumps are interpreted as being the result of heavily over- 
lapped pulses in §5.2. Finally, in §pl we discuss how our 
results impose important constraints on the current phys- 
ical models and in particular how they will be relevant in 
comparing external versus internal shock models. 

2. THE RELATION BETWEEN HARDNESS AND INTENSITY 

The relation between the hardness and the intensity, 
during the active 7-ray phase of a GRB, has been well 
investigated. It has been shown that there is no ubiqui- 
tous trend of spectral evolution that can characterize all 
bursts; several types of behavior exist. Norris et al. (1986) 
found that the most common trend of spectral evolution 
is a hard-to-soft behavior over a pulse, with the hardness 
decreasing monotonically as the flux rises and falls. This 
conclusion was also arrived at by Kargatis et al. (1994, 



hereafter K94). A few cases exhibited soft-to-hard and 
even soft-to-hard-to-soft evolution. Band (1997) studied 
the data from the four, high time-resolution channels from 
the Large Area Detectors (LADs) of BATSE. The spectral 
evolution was analyzed by auto- and cross-correlating light 
curves from the different energy channels. Most bursts in 
the sample showed a hard-to-soft behavior. 

There are also bursts that do not seem to exhibit any 
HIC at all, with an apparently chaotic behavior. The 
main conclusion drawn by Laros et al. (1985) and Jour- 
dain (1990), was that there did not exist a HIC between 
the spectral evolution and the light curve in their sam- 
ples (using PVO and APEX data, respectively). Over the 
whole GRB, there often does not exist any pure corre- 
lation, even though the tracks in the hardness-intensity 
plane are confined to an area from hard and bright to soft 
and weak, indicating an overall trend of increasing lumi- 
nosity with hardness (K94). A seemingly chaotic behavior 
in that plane may be the result of a superposition of several 
short hard-to-soft pulses that cannot be resolved. Various 
types of trends have also been seen in a single GRB (e.g., 
Hurley et al. 1992). The variety of behaviors are also de- 
scribed in Band et al. (1993) and Ford et al. (1995). 

Furthermore, there is another behavior in which the in- 
tensity and the hardness track each other. This behav- 
ior is less common than the hard-to-soft trend and was 
first noted by G83, who described it quantitatively. They 
found a power law relation between the instantaneous en- 
ergy flux, F, and the energy parameter Eq derived from 
modeling the photon spectra using Ne(E) — E a e~ E ' E °, 
which serves as a measure of the hardness, i.e., 
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The power law index 7 was found to have a typical value 
of 1.5-1.7. This value is sometimes referred to in the lit- 
erature as the correlation index. 

This analysis was criticized by several workers, includ- 
ing Laros et al. (1985), Norris et al. (1986), and K94. It 
was speculated that the correlation could possibly be an 
artifact of the way the hardness was derived from the two- 
spectral-channel count rates. Furthermore, G83 excluded 
the hard initial phase of the bursts. Ford et al. (1995) 
suggested that the low time-resolution may result in the 
initial, hard behavior being missed. On the other hand, 
K94 confirmed the description of the HIC made by G83 
(eq. pj), i.e., a power law model of the hardness-intensity 
correlation (hereafter denoted by PLHIC). They found a 
power law correlation in approximately half of their cases, 
but with a substantially wider spread, 7 = 2.2 ± 1.0. Fur- 
thermore, Strohmayer et al. (1998) investigated the evolu- 
tion of the peak energy versus the energy flux in the Ginga 
data and found the PLHIC to be valid here too, with, for 
instance, 7 - 3 for GRB 890929 (in the 2-400 keV energy 
range). 

The present paper concerns mainly the power-law, 
hardness-intensity correlation during the decay phase of 
pulses, which is a common behavior. Kargatis et al. (1995) 
found the PLHIC in 28 pulse decays from 15 bursts, out 
of a total of 26 GRBs with prominent pulses studied. The 
distribution of the correlation index peaks at 1.7 and has 
a substantial spread. A large spread in the PLHIC index 
was also found by Bhat et al. (1994), who studied 19 time 
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structures with short rise times (< 4 s) and slow decays, 
and concluded that most had a good correlation between 
the hardness and the intensity. The value of 7 varied from 
1.4 to 3.4. 

Ryde & Svensson (2000a, 2000b) derived the conse- 
quences of combining the power law model of the HIC 
and the exponential model of the Hardness-Fluence Cor- 
relation, for the decay phase of GRB pulses. They found a 
self-consistent, quantitative, and compact description for 
the temporal evolution of the pulse decay phase. It was 
shown that, assuming the adopted models are valid, the 
total photon flux must be oc 1/(1 + t/r), where the time 
t is taken from the start of the decay and r is a time con- 
stant that can be expressed in terms of the parameters of 
the two empirical correlations. 

3. DATA AND METHODS 

3.1. Observations 

This work is based on the data taken by BATSE on 
board the CGRO (Fishman et al. 1989). It consists of 
eight modules placed on each corner of the satellite, giving 
full sky coverage. Each module has two types of detectors: 
the Large Area Detector (LAD) and the Spectroscopy De- 
tector (SD). The former has a larger collecting area and is 
suited for spectral continuum studies, while the latter was 
designed for the search of spectral features (lines). For 
our spectral analysis we used the high energy resolution 
(HER) background and burst data types from the LADs 
having 128 energy channels. The burst data have a time- 
resolution of multiples of 64 ms. The CGRO Science Sup- 
port Center (GROSSC) at Goddard Space Flight Center 
(GSFC) provides these data as processed, high-level prod- 
ucts in its public archive. The data are available for all 
the detectors that triggered on the bursts (often 3 or 4 of 
the detectors closest to the line-of-sight to the burst loca- 
tion). Models of the relevant Detector Response Matrix 
(DRM) for each observation are also provided (Pendleton 
et al. 1995). The eight modules of BATSE allow the local- 
ization of the GRB, needed for the determination of the 
DRM, since it is dependent on the source-to-detector axis 
angle. Finally, we used, for visual inspection of the light 
curves, the so called concatenated 64-ms burst data, pro- 
vided by GROSSC, which is a concatenation of the three 
BATSE data types DISCLA, PREB, and DISCSC. 

3.2. Sample Selection 

To construct a complete sample of strong bursts, we 
started by selecting the bursts in the Current BATSE 
Catalog 1 , up to GRB 990126 (BATSE trigger number 
7353), for which it is possible to measure peak fluxes. 
These are approximately 80% of the totally 2302 ob- 
served. Data gaps and/or missing data types are the 
reason why the peak flux for some bursts are not found. 
The Current BATSE Catalog is preliminary, but bursts 
up to GRB 960825 (trigger 5586) are published in the 4th 
BATSE Catalog (Paciesas et al. 1999). The threshold we 
chose for accepting a burst was set to a peak flux (50-300 
keV in 1.024 s time resolution) of 2 photons s" 1 cm~ 2 . 
This selection resulted in a set of 420 bursts. 

This set was examined visually, case by case, using the 
concatenated 64 ms time resolution data. We searched for 

1 http://gammaray.msfc. nasa.gov/batse/ 



bursts containing long pulse structures with a general "fast 
rise-slow decay" , often referred to as "fast rise-exponential 
decay" (FREDs). No analytical function describing the 
pulse shape was assumed. The reason for using such a 
loose definition is to have a sample that is independent 
of any preconceived idea of the pulse shape. Examples of 
pulses that are subsequently included in this broad sample 
are shown in Figure |lj. 

The identification of the pulse structures may, in some 
cases, be disputable. A few algorithms for identifica- 
tion have been introduced by, e.g., Li & Fenimore (1996), 
Pendleton et al. (1997), and Scargle (1998). Norris et al. 
(1996) developed a method to identify pulses based on as- 
suming stretched exponential pulse shapes. This method 
was used by, e.g., Crider et al. (1999) to separate pulses, 
some of which were heavily overlapped. However, all these 
methods depend on various assumptions. Our sample suf- 
fers, on the other hand, from some subjectivity, as we se- 
lect our pulses visually. However, this should not affect the 
results, as these are not strongly dependent on the details 
in this selection process. 

For the time-resolved spectroscopy, we will use a high 
signal-to-noise ratio (S/N), which will lead to light curves 
consisting of only a few broad time bins. We need as many 
time bins as possible to study the HIC and to arrive at re- 
liable results for the PLHIC index, as we will be determin- 



ing its distribution (sec §4.2 for details). For this purpose 
we adopt the criterion that the decay phase of the pulses 
should have at least 5 time bins with S/N = 30 to be 
included in the study. 

It should be noted that the first and the last time bins 
could partially cover a time interval that is outside the ac- 
tual HIC behavior of the decay phase. For the pulse decay, 
the first time bin studied corresponds to the peak of the 
pulse, and therefore this time interval might include the 
transition between the rise and the decay phase. The last 
time bin might also be a transition interval, i.e., the valley 
before the rise phase of the next pulse. In the 5-time-bin 
cases, bearing these qualitative arguments in mind, there 
are, in the worst case, only 3 central data points that are 
unaffected and thus are more certain to correctly define 
the PLHIC index. 

The visual inspection resulted in a set of 66 bursts, that 
is ss 3% of the original BATSE catalog. These bursts are 
presented in Table 1 and they constitute our main sample, 
which we have striven to make as unbiased and extensive 
as possible. In the Table, the bursts are denoted by both 
their BATSE catalog and trigger numbers. The detector 
from which the data were taken and the number of time 
bins (nbins) selected for fitting are also presented. The 
pulses studied are identified by the time t max when their 
count rate is maximal (t — is the trigger time). 

3.3. Spectral Modeling 

The central part of the analysis was performed with the 
WINGSPAN package, version 4.4.1 (Preece et al. 1996a). 
The spectral fitting was done using the MFIT package, 
version 4.6, running under WINGSPAN. We always chose 
the data taken with the detector which was closest to the 
line-of-sight to the location of the GRB, as it has the 
strongest signal (see Table fOfor the individual cases). The 
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broadest energy band with useful data was selected, often 
25-1900 keV. A background estimate was made using the 
HER data, which consist of low (16-500 s) time resolution 
measurements that are stored between triggers. The light 
curve of the background, during the outburst, was mod- 
eled by interpolating these data, roughly 1000 s before and 
after the trigger, with a second or third order polynomial 
fit. 

To perform detailed time-resolved spectroscopy it has 
been shown that a S/N ~ 45 is needed (Preece et al. 
1998). The aim of our spectral analysis, of every time bin, 
is mainly to determine the peak energy, E p k, as a mea- 
sure of the hardness, and allowing a deconvolution of the 
count spectrum to find the energy spectrum. Therefore, 
we accepted a lower S/N, sometimes down to 30, but we 
checked that the results are consistent with higher S/N ra- 
tios. This gives us the possibility to study the burst pulses 
with higher time-resolution, which is of great importance 
for our study. 

The background-subtracted photon spectrum, Ne{E), 
for each time bin, is then determined using a forward- 
folding technique. An empirical spectral model is folded 
through the model of the Detector Response Matrix and 
is then fitted by minimizing the % 2 (using the Levenberg- 
Marquardt algorithm) between the model count spectrum 
and the observed count spectrum, giving the best-fit spec- 
tral parameters and the normalization. The spectra were 
modeled with the empirical function (Band et al. 1993): 



( , ( AE a e~ E / E ° ii(a-p)E >E 



(2) 



where E is the energy, Eq is the e— folding energy, a and 
/3 are the asymptotic power law indices, A the amplitude, 
and A' has been chosen to make the photon spectrum 
Ne {E) a continuous and a continuously differentiable func- 
tion through the condition 



A' = A [{a - (3)E ] a 
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The power law indices were always left free to vary if noth- 
ing else is stated. The peak energy, E p ^, at which the 
-E-FE-spectrum (Fe = ENe) is at its maximum, was used 
as a measure of the spectral hardness instead of Eq . They 
are related by -E p k = (2 + a)Eo and a peak exists only 
when (3 < — 2 < a. The fitting procedure has then 4 free 
parameters: A, a and /?, and _E p k- The photon spectrum 
arrived at is model-dependent. However, as equation Q2J) 
often gives a good model of the spectra, the photon spec- 
trum found by deconvolving the count spectrum should 
correspond well with the true photon spectrum. For every 
time bin, the instantaneous integrated energy flux F was 
found by integrating the modeled energy spectrum over 
the available energy band of the detector. 

These procedures have now given us, beside the spec- 
tral parameters, a data set of peak energies and energy 
flux values for every time bin of the pulse decays. Note, 
however, that due to the low energy cut-off of the detec- 
tor, in practice only E v ^ values > 40 keV can be reliably 
determined. For this reason, it is often the case that a few 
time bins, at the end of the pulse decays, are not used in 
our studies of the hardness-intensity temporal evolution. 



3.4. HIC Analysis Method 
3.4.1. Statistical Analysis 

Different statistical approaches can be used to test 
whether a power law model is adequate to describe the 
hardness-intensity data, and to determine its best-fit pa- 
rameters and their likely uncertainties. Here, we once 
again used the minimization of the merit function x 2 ■ If 
the uncertainties of the measurements are known, a model- 
independent determination of the goodness-of-fit can be 
obtained. This is done typically by calculating the proba- 
bility Q of exceeding by chance the \ 2 minimum obtained 
(see, e.g., Press et al., 1992). Data points were weighted 
using the uncertainties in the -E p k(£) estimation, which are 
usually much larger than the intensity uncertainties. We 
found in too many cases unrealistic values of Q ~ 1. The 
most likely explanation is that the uncertainties derived 
from the spectral fitting are overestimated. Their use- 
fulness is therefore doubtful. One possible cause of this 
could be the strong correlation observed between some of 
the Band et al. function parameters, particularly between 
a and A (see, e.g., Lloyd & Petrosian 2000). Such a cor- 
relation could partly be due to an artifact of the fitting. 
However, to distinguish this from a physical relation be- 
tween the parameters, a more detailed study has to be 
undertaken. 

The power law model can be linearized by taking the 
logarithms of the intensity and the hardness measure. The 
power law index then becomes the slope of the linear re- 
lation, and ordinary linear regression can be used to fit 
the data. After this transformation, in most cases, the 
-E p k-measurement uncertainties are approximately equal 
and, since the detector energy channels are approximately 
equally spaced on a logarithmic scale, more symmetric er- 
ror bars are expected. In addition, the HIC data points 
are more evenly distributed. Thus, assuming that the un- 
certainties are exactly equal, i.e., using the unweighted 
method, we can obtain the coefficient of determination, 
R 2 , i.e., the square of the linear correlation coefficient. 
For any given R 2 , Pn(R 2 ) is the probability that N mea- 
surements of two uncorrelated variables would give a co- 
efficient larger than R 2 . This probability is a measure of 
how significant the linear correlation is. 

The tests R 2 , x 2 , and the like, do not have information 
about the temporal order of the data points, so obviously 
one should not infer anything about the temporal evolu- 
tion of the HIC from these statistics. However, most pulse 
decays in our sample that have a very good PLHIC also 
show a good tracking behavior, i.e., the temporal evolution 
in the hardness-intensity plane is monotonic, aside from 
random fluctuations that can be attributed to the mea- 
surement uncertainties. Nevertheless, this generalization 
must be taken with caution in cases with low correlation 
coefficients. Another important aspect not measured by 
these tests is whether the residuals show any particular 
trend or feature. In this respect, we analyzed in detail the 
cases with many time bins and concluded that the data 
are consistent with the power law model. 

To fit the data directly one can employ a non- 
linear regression numerical algorithm (with the Levenberg- 
Marquardt method). The outcome is essentially the same 
as in the linear regression and will not be presented here. 

In conclusion, the difference between the indices ob- 
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tained using weighted and non-weighted fittings was, in 
most cases, within the estimated uncertainties. No sig- 
nificant difference was obtained either from our various 
statistical analyses of the two data sets. We will show the 
results derived from the non- weighted, linear fittings, but 
we checked that our conclusions are independent of this 
choice. 

3.4.2. The (p-method for studying the HIC 

The limited spectral coverage of the detector used might 
affect the assigned measure of the bolometric flux, i.e., the 
energy-integrated flux, as a substantial fraction of the flux 
could be lost, especially when the spectrum has a broad 
shape or peaks close to the boundaries of the detector. A 
second problem, when one aims at studying the correla- 
tion within single pulses in the light curve, arises from the 
fact that the observed spectra may contain contributions 
from other pulses, for instance, previous pulses which still 
contribute with soft photons, and unresolved, overlapping 
pulses. Furthermore, additional, separate soft components 
(Preece et al. 1996b) could also alter the measured flux 
value and thus weaken the correlations. All this can affect 
the analysis by changing the shape of the spectrum. 

This motivated us to introduce a new representation of 
the HIC, which might resolve some of these complications 
(see also Ryde et al. 2000). The value of EF E at P pk can 
be used as a representation of the energy- integrated flux as 
it, under some circumstances, is proportional to the total 
flux. This quantity will be denoted by ip (see Fig. Ha) and 
the PLHIC can be studied as 



if oc EP 



pk' 



(4) 



where r\ is the new PLHIC index. This discussion is lim- 
ited to the cases where the EFe peak actually exists within 
the detector band [E min , E max ], which most often is the 
case (Band et al. 1993). In the most common case, where 
-Emax > (a — /3)E p k/(a + 2), the proportionality between 
cp and F — J F E dE becomes 



F e («+2) 

= K(a, p, J/minj 2/max) — 7 nlrv-l-9 X 

ip (a + 2) a+z 

[T(a + 2){P{a + 2, a - 0) - P(a + 2, y min )} 



(5) 
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where y min = (a + 2)E min /E pk , and y max = (a + 
2)E max / 'Epk. r(a) and P(a,y) are the gamma func- 
tion and the incomplete gamma function, respectively 
(see, e.g., Press et al. 1992). In the case that a and 
depend weakly on -E" p k, the only £" p k dependence in 
k(a, 0, y m i„, y max ) is in y min and y max . In particular, when 
the flux integration is chosen to be over the whole energy 
range from to oo, there will be no dependence at all. 
Therefore, ip should be a better representation of the bolo- 
metric flux for the study of the HIC. 

4. RESULTS OF THE ANALYSIS 

4.1. Correspondence between the F and ip 

Before discussing the results of the HIC analysis, we 
will compare the HIC relation as given by equation (|J), 



bu t us ing the parameter E pk as a measure of the hardness 

(§© 

F*E; k , (6) 

and the new representation given by equation (0). For 
this comparison, we use a subset of 47 pulse decays (in 
39 bursts) from the sample, which have been selected as 
having good PLHICs, i.e., the cases having a relative un- 
certainty An/r] < 0.15 and marked with an A in Table M 
(see § |4.2| for a discussion of this choice). For each pulse 
decay we analyzed the spectral evolution and determined 
both the 7 and the r\ values. When comparing the coeffi- 
cients of determination given by the two methods we found 
that in 83% of the cases Rl is greater than Rl. To find 
whether the observed difference is significant, we made use 
of Fisher's z-transformation 



Z= V— 

2 \\-R 



(7) 



where each z is approximately normally distributed, with 
a standard deviation a(z) ss l/\/ N — ~3, and N bein g the 



number of data points (see, e.g., |Press et al. 1992J ). In- 



dividual differences are not significant, but the difference 



between the respective mean values, 



was found 



to be positive at a significance level of P-value = 0.015 
(the P-value is the probability that the value of the test 
statistic is as extreme as it is, with the null hypothesis be- 
ing true). This implies that the new ip-method does give 
better correlations. 

The relation between the two PLHIC indices 7 and 77 for 
all pulse decays in the subset A is shown in Figure g. A lin- 
ear fit to these data gives 7 = (0.96±0.09)r7+(0.08±0.16). 
This shows that there is a good average correspondence be- 
tween the two methods. All these results convinced us to 
use the </j-method in the subsequent studies. 

4.2. The PLHIC index of Single Pulses 

The results of the fitting of the bursts in our main 
sample are presented in Table @. Apart from the power 
law indices 77 and 7 obtained, it shows B% and its associ- 
ated probability Pn(B1), R'i, and the relative uncertainty 
At?/?/. 

From these bursts we select and study the subset A for 
which the hardness-intensity relation is well fitted by a 
power law model. This decision is not trivial and the con- 
cern was to set a reliable rejection level to select those 
cases that are consistent with the model. The proba bility 
Pn(R 2 ) could be used for this purpose (see §3.4.1). In 
practice, the relative uncertainty of the slope in the linear 
regression, either A77/7/ or A7/7, gives a more convenient 
measure (at least when the range of slope values is not 
close to zero). We found empirically that, for our set, a 
rejection level of P/v(P^) < 0.001 (usually considered as 
highly significant) is approximately equivalent to selecting 
cases with Arj/rj < 0.15. This level was chosen to define 
the subset A in Table 0. We used the r\ indices in our 
select ion because the </?-method gives better correlations 
(see §|4.l| ), but almost the same subset would be reached 
using A7/7 < 0.2. The selection of the rejection level was 
a trade-off between choosing a high level and reducing sig- 
nificantly the sample size, or a low one that would increase 
the statistics but may allow cases that are likely also to 
have large systematic errors. 
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4.2.1. The Distribution of PLHIC Indices among Bursts 

To study the distribution of some parameter associated 
with GRB pulses, one should consider that pulses within 
the same burst may not be independent. Furthermore, 
they may have a different distribution of the parameter 
values considered as compared to their distribution for 
pulses over the full burst sample, i.e. the general distribu- 
tion. We will show below that at least the latter is true 
for the 7 and r\ indices. In that case, a small sample can 
be easily biased by taking many values from a multi-pulse 
burst. Therefore, to estimate the general distribution of 
the PLHIC indices when two or more pulses were mea- 
sured in the same burst, only one was taken. To select it, 
we chose the one that shows the best correlation in terms of 
the probability Ppf(Bl). This method is consistent with 
the fact that very often only one pulse, among many in 
the burst, is found suitable for this analysis, i.e., we are 
already disregarding noisy and therefore poorly correlated 
cases. On the other hand, since indices are very similar 
within bursts, as will be found below, other criteria, such 
as a random selection, produced no significant differences. 

We studied how the resulting distributions depend on 
the chosen rejection level Ary/77. When trying to estimate 
the underlying probability distribution of any sample, care 
must be taken that the measurement uncertainties do not 
contribute significantly. Allowing cases with larger rela- 
tive uncertainties, i.e., worse PLHIC cases, results in al- 
most identical mean values but larger standard deviations. 
This is to be expected, since the data have now a larger 
intrinsic dispersion that is added to (convolved with) the 
real distribution. For later use in our studies, it is impor- 
tant to set a fairly high rejection level, in order to have a 
good estimate of the real dispersion. 

Thus, taking from the main sample one pulse per burst 
and only the cases with A 77/77 < 0.15, defines the subset 
of 39 pulses marked with B in Table @, for which we will 
study the general distribution of the indices (note that 
B C A) . Using the F- method, the mean value of the PL- 
HIC index is 7 = 1.9 ±0.1 with a standard deviation of the 
sample er 7 = 0.72 ±0.08. We also provide the expected de- 
viations of these estimators to facilitate the comparisons. 
Using the 93-method the mean value and dispersion for the 
PLHIC indices are fj = 2.0 ±0.1 a nd a„ = 0.68 ±0.08. The 
weighted fittings mentioned in (§ 3.4.1) produce very sim- 
ilar results. For example, the values obtained with that 
method are fj w = 2.15 ± 0.10 and cr r , w = 0.80 ± 0.09. The 
differences are of the same order as the estimated uncer- 
tainties. 

These results should only be compared with the ones 
obtained in studies of the decay part of pulses. As men- 
tioned above, Kargatis et al. (1995) found a central value 
of 7 ~ 1.7, which is somewhat lower than ours. However, 
their sample of 28 pulses was taken from 15 multi-pulse 
GRBs. This might produce a slower convergence to the 
real mean value. In fact, if we allow more than one pulse 
per burst in our statistics, we get a 7 = 1.8 ± 0.1 for 47 
pulses from 39 bursts (subset A in Table g). Note that 
although in many bursts only one pulse was found suit- 
able for our studies, it does not mean that there were no 
other pulses present. A real distinction between single and 
multi-pulse bursts would be difficcult to achieve, especially 
as we do not know what a single pulse really looks like. 



In Figure g, histograms of the subset B of measured 
PLHIC indices, 7 and 77, are displayed. The sample is too 
small for a detailed study of the underlying distributions. 
However, for later use in this work (§ |Q| ), we test the 
assumption of a normal distribution for both indices. We 
use the Geary test of normality (see, for instance, Devore 
1982), a simpler but more powerful test procedure than 
the x 2 test, that is specific for this purpose. It is based on 
the ratio of the average absolute deviation to the square 
root of the average square deviation, given by 



U = 



X\/n 



T,i(Xi-X) 2 /n 



(8) 



For samples larger than 20 the Geary test can be approx- 
imated by a two-tailed Z-tcst, where the null hypothesis 
is that the underlying distribution is normal and the stan- 
dardized U is 

For the r\ and 7 samples, high P- values are found, 0.73 
and 0.34 respectively, implying that the approximation of 
a normal distribution can be reliably used. 

4.2.2. The Excluded Cases 

We found, using the R 2 statistics, that 57% of the main 
sample of pulse decays (subset A) could be described by a 
power law at a highly significant confidence level. At the 
level of significance of only Pm{R 2 ) < 0.01 (approximately 
equivalent results would be obtained using A77/77 < 0.2) 
the number increases to about 75%. Since the R 2 coeffi- 
cient is not a very sensitive test, these percentages should 
be taken as approximative. 

There are several possible causes for some of the cases 
showing poor power law correlations. We could often rec- 
ognize problems that were affecting our analysis of the 
HIC in a systematic way, making uncertain any conclu- 
sion about the applicability of the power law, or any other 
model for that matter. In some cases, for instance, the 
-Ei^-spectrum peaks at higher energies than the observed 
band. However, the numerical algorithm finds an E v ^ 
value within the energy window using > — 2. As the 
whole spectrum evolves towards lower energies, at some 
point the absolute maximum enters the observation win- 
dow but the overall evolution of the HIC is poorly fitted 
by a power law. 

Although the Band et al. function proved to be a very 
flexible model to describe most GRB spectra, it presents 
difficulties when trying to fit some particularly broad 
shaped cases (see, e.g., Ryde 1999b). In a graphical EFe 
vs. E representation, a typical zig-zag evolution is ob- 
tained, accompanied by opposite changes in a and j3. A 
similar effect occurs when either of the latter two parame- 
ters is loosely constrained by the data. Although in these 
cases an overall track behavior is apparent, PLHIC fits re- 
sult in poor R 2 coefficients, and residuals show an uneven 
dispersion. 

Furthermore, the selected pulse decay could also turn 
out not to be a single pulse decay after all, even if it ap- 
pears to be. Overlapping, "hidden" secondary pulses in 
the light curve produce a characteristic effect in the HIC 
evolution that, under some circumstances, can be recog- 
nized and used as a pulse identification aid. This will be 
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discussed in detail in § [4.3.1 . Finally, some cases could 
actually have another HIC behavior. An investigation of 
alternative models will be presented elsewhere (Ryde & 
Svensson 2000c, in preparation). 

4.3. The PLHIC Index of Multi-pulse Bursts 

We will now study the distribution of the PLHIC indices 
within a single burst. In order to increase the sample, we 
look for additional pulses, with 4 time bins, in the GRBs 
listed in Table pi. The expanded sample of multi-pulse 
bursts includes 4 of these shorter pulses, but in each case 
there is a longer pulse to compare with. The pulses from 
GRB 950104 (t rigger 3345) are also included, based on the 
discussion in § 4.3.1. In this way, we increase the number 
of multi-pulse cases from 11 to 15 GRBs. They all show 
at least moderately good correlations. 

Apart from these cases, a set of pulses, apparently sin- 
gle, but having a track jump in the HIC will be studied in 
detail. They are noted with TJ in Tabled and can be inter- 
preted as consisting of several pulses (see also Borgonovo 
& Ryde 2000). 

4.3.1. Bursts with Several Well-separated Pulses 

The sample of GRBs which have two or more well- 
separated pulses is shown in Table 0. As an example, 
we present the case of GRB 921207 (trigger 2083) in Fig- 
ure m This burst ha s been exam ined by man y workers 
(e.g-T IFord et al. 1995| , Crider et al. [19984 11999| , and Ryde 
& Svensson 1999a ]and is one of the best BATSE cases 
available for this type of study, with two bright and long 
pulses. In Figures He and d, the hardness-intensity evolu- 
tion is shown. The decay phases of both pulses show very 
good correlations and the power law indices are equal to 
within the estimated uncertainties. ^From the point of 
view of the PLHIC, the first pulse structure is consistent 
with being a single pulse. This is also the conclusion of 
Ryde & Svensson (1999) who studied various aspects of 
the spectral/temporal evolution of GRB pulses. However, 
with the pulse identification method proposed by Norris 
et al. (1996) this pulse structure can be modeled as a su- 
perposition of two individual, overlapping, stretched ex- 
ponential pulses (see Crider et al. 1999). In the present 
spectral analysis there is no indication of this, although 
one cannot rule it out. 

As another example, the overall evolution of 
GRB 970420 (trigger 6198) is shown in Figure |. While 
the spectral evolution during the rise phases does not seem 
to follow any clear trend, again in this case we observe 
a good PLHIC during the decay phases. The slopes are 
equal to within the estimated uncertainties (see Table 0) . 

In Figure § the light curve of GRB 950104 (trigger 3345) 
is shown. The section after the peak at t ~ 6 s can be in- 
terpreted as having two pulses which overlap closely. An 
indication of the existence of a secondary pulse can be 
found as follows. The assumption has to be made that 
the rise and decay phases of individual pulses are strictly 
monotonic and then we calculate whether the deviation 
can be attributed to the Poisson noise in the signal. For 
that purpose, a simple test would be to find the signifi- 
cance of the count difference between the suspected peak 
and the nearby valleys (see Li & Fenimore 1996): 



where N a is the number of a standard deviations. In the 
case of GRB 950104, we found a significance N a = 3.2, 
within the interval 3 < N a < 5 recommended for the test, 
i.e., marginally significant. Figures a and b show two 
count time history plots of this burst; one with a time 
binning with S/N = 40 and the other with the 64-ms 
time resolution of the concatenated data. The spectral 
evolution of the peak energy corresponding to the time 
binning in panel (a) is also shown in panel (c). Initially, 
the second decay has approximately the same PLHIC in- 
dex as the first one. It then makes a change into a par- 
allel track. Checking the 64 ms count rate time history, 
one can see that this track jump coincides with the pres- 
ence of the small secondary peak (marked with an arrow) 
which is hidden in the low time resolution binning. Even 
though the significance of the peak is somewhat marginal, 
the agreement with the slope of the first pulse track makes 
the case more compelling for the presence of overlapping 
pulses. The power law indices obtained for the tracks are 
?/i = 0.77 ± 0.05, 772 = 0.79 ± 0.03, and ?y 3 = 0.76 ± 0.05, 
equal to within the estimated uncertainties. 

Note that the -Enk uncertainties derived from the spec- 
tral fittings (Fig. gc) are heavily overestimated. This is a 
common charact eris tic of the track jump cases, and it will 



be discussed in §5.1 



4.3.2. Apparently Single Pulses 

The three cases presented here seem to be single pulses 
as seen from the 64 ms light curves. These pulses are 
taken from GRB 910927 (trigger 829), GRB 960912 (trig- 
ger 5601), and GRB 970925 (trigger 6397), respectively 
and they are shown in Figures [| and [7]. All three behave 
similarly to the case GRB 950104 discussed above. We 
have modeled their decay evolution with two power laws 
and a track jump. The result of the fittings are shown in 
Table |. 

Although there is no clear indication of secondary peaks, 
in the case of GRB 970925, a small "bump" in the light 
curve (more easily noticed in a log-linear graph) coin- 
cides with the track jump. The measured indices, ?/i = 
0.77 ±0.15 and 772 = 0.86 ±0.05, are again equal to within 
the uncertainties, as in the other two cases. 

In the case of the pulse decay in GRB 910927, and only 
in this case, the fitted results are obtained with the pa- 
rameter j3 frozen for all points. Although the qualitative 
behavior is the same, it produces a sharper transition be- 
tween the tracks. Note that the spectral evolution (Fig. 0c) 
during the decay of the first pulse, given by the bins 2-4, 
seems to follow the same track as do the points numbered 
14-16. It could therefore be that the last points, in fact, 
belong to the decay of the first pulse and that the pulse 
peaking at bin 9 has a fast decline and is superimposed in 
the middle of the dominating, long first pulse. 



4.4. The Distribution of the PLHIC Index within GRBs 

An interesting question is how the distribution of the 
HIC index within a bur st c ompares to the general distri- 
bution as presented in §4.2. Our set of multi-pulse bursts 
was introduced in §4.2 and listed in Table g. It consists 
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bles 



a pulses taken from lo UHBs. We can see, trom la- 
g and |[ that the PLHIC indices are, in most cases, 
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very similar. To estimate the probability of having similar 
indices within bursts purely by chance, it will be assumed 
first that the set of random variables Xi £ {77} (or {7}) 
have approximately a normal d istrib ution Xi ~ N(fj,a v ), 



as has been justified in section §4.2.1, Let Y = X. 



■Xjbe 

the difference between two indices from the same burst. 
If they are independent, then Y ~ N(0,^/2a v ). A hy- 
pothesis test can be performed to see whether cr 2 , = 2er 2 
is the true variance. Thus, choosing the null hypothesis 
Hq : cr 2 = Cq , we will use first the test statistic 



2= (n-l)s 2 



(11) 



where s is the standard deviation and n — 1 is the number 
of degrees of freedom. 

This test has the advantage of using a standard proba- 
bility distribution function, but special care must be taken 
when considering the cases where 3 (or more) pulses have 
been measured within a burst, let us say A, B, and C. We 
can form three differences, but obviously one of them will 
not be independent of the others. An approximation could 
be done adjusting the number of degrees of freedom, but 
then equation (O) will not be strictly valid. Numerical 
simulations convinced us that, at least in this case, this 
is not a good approximation. The results shown in Ta- 
ble p| from the x a test were obtained selecting pairs (A, B) 
ana (B, C). We verified that our conclusions are not af- 
fected by this choice. The dispersions expected assuming 
independency are y/2a 7 — 1.02 and y/2a v = 0.96, and the 
standard deviations obtained are s 7 = 0.39 and s n = 0.30, 
respectively. In all cases, the P-values are very small for 
both indices 7 and r). Therefore we can confidently reject 
Hq in favor of the alternative hypothesis, i.e., the distribu- 
tion of PLHIC indices over a single burst is narrower than 
the one observed over the whole sample of GRBs. 

In Figures @a and b we illustrate how all possible index 
differences, taking pulses in temporal order, compare with 
the normal distribution assumed in Hq. We had to run a 
numerical simulation to find the non-standard probability 
distribution for the variance of all index differences. We 
define the test statistic 
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where N is the number of bursts in the sample and Mi is 
the number of measured pulses within burst i. This test 
is obviously independent of the order in which the differ- 
ences are taken. Using a Monte Carlo method, we found 
the probability density function of T a assuming, as before, 
that each index Xij follows the general distributions found 



in § 4.2.1, and using the normal approximation. The re- 
sults of this second test are presented in Table [| Again, 
P-values are very small for both indices, confirming the 
conclusions of the previous test. The purpose of using first 
the xi l es t as an approximation is that, having a known 
analytical solution, it can estimate very low probabilities 
without being limited by computational time. 

We have now established that the distribution of indices 
among pulses within a burst is narrower than the one be- 
tween pulses of different bursts. Figures Ha and b show 
that most uncertainties in our sample are much smaller 



than the dispersion of the general distribution. Note also 
that most error bars cover the origin, suggesting that in- 
dividual index differences are not significant. We will now 
consider whether the index dispersion within a burst could 
be attributed almost entirely to the stochastic uncertain- 
ties. Figure gjc and d show the standardized residuals of 
all possible differences Y t between pulses of the same burst 
in our sample. The proportion of data points within inte- 
ger multiples of a is consistent with a Gaussian distribu- 
tion, except for an "outlier" in the 77 set, at number 14. 
This point and point number 12 are both differences taken 
against the second pulse of GRB 950624 (pulse 3648b; see 
Table ph which is a fairly weak but long pulse. It precedes 
and slightly overlaps a strong pulse (3648c), which could 
produce a systematic error in our measurements. Note 
also, however, in Table || that comparing the estimated 
index uncertainties almost always 2 ;> A7/A77 ;> 1. One 
exception is this case where A7/A77 w 8. This is why for 
the 7 differences (Fig. ||c) the corresponding point is within 
±lcr. The cause of the large deviation is most likely due 
to a casual alignment of data points that resulted in an 
underestimation of the A77 uncertainty. This fact stresses 
the importance of having high selection requirements. If 
we had tried to expand our sample to include even weaker 
or shorter pulses, with fewer time-bins, it would have been 
at the cost of increasing the risk of systematic errors. Since 
these are difficult to identify and quantify, they can easily 
distort the statistics. 

Assuming that the indices are constant within a burst, 
we will estimate the probability to exceed the observed dis- 
persion. The uncertainties will be assumed to be drawn 
from a normal distribution. Therefore the index differ- 
ences Y £ {Xij — Xik} should follow Yi ~ N(0, <Xj), where 
each of is found by adding the corresponding Xi variances. 
We define, in a similar way as before, a test statistic 
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where the summation considers all possible differences. 
Since not all differences are independent, we have to solve 
it numerically. On the other hand, if we restrict the sum 
to independent terms 



xNE E 1 ^ 



Xij+l) 



;-i ,-1 a ij +a ij+l 



(14) 



the test becomes a standard x 2 statistic. 

The results of these last two tests are summarized in the 
second part of Table [|. For the 7 indices we found rela- 
tively low P-values for both tests, but we consider that, 
at the 0.07 probability given by the more precise test (T ), 
the hypothesis of invariance of the 7 index is marginally 
acceptable. The results concerning the 77 indices depend 
on the inclusion or rejection of pulse number 3648b men- 
tioned above. The outlier term completely dominates the 
sums in equations ( |14| ) and (O). Its inclusion gives P- 
values <C 10~ 5 for both tests. Obviously, a 9<r difference 
can not be explained in terms of Gaussian uncertainties. 
On the other hand, rejecting this data point, we obtained 
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high P- values, implying that the observed differences be- 
tween 77 indices are most likely due to the measurement 
uncertainties. 

Finally, considering again Figure @, it is apparent that 
while the 77 differences have approximately a zero mean, 
the 7 differences tend to be negative, i.e., later pulses have 
generally larger indices. Taking the averages of the stan- 
dardized residuals, the means are —0.037 and —0.49 for 
the 77 and 7 differences, respectively. We calculated, in the 
7 case, whether the non-zero mean deviation is significant. 
For that purpose a similar statistic to equation ( |l3| ) was 
used, but without squaring the terms. We found a 0.03 
probability (two-sided level) to exceed the observed value 
assuming a zero mean, indicating a significant deviation. 

5. DISCUSSION 
5.1. Comparison between the F and (p methods 

The introdu ction of the (^-method for HIC studies was 
motivated in § 3.4.2 . It was shown (Fig. ||) that there is a 
good average correspondence between the indices 7 and 77. 
The scatter around an exact correspondence was to be ex- 
pected. As mentioned above, additional flux components 
can be contributing to the estimated flux affecting the re- 
sulting PLHIC index. The ^-method is less dependent 
on this. Furthermore, the exact equivalence between the 
power law indices is valid only in the case that a and do 
not vary with -E p k, which they do in many bursts. Further- 
more, the flux F is found by integrating the deconvolved 
spectrum and the deconvolution is model-dependent. This 
could introduce additional scatter into the PLHIC relation. 

We have used both F and <p in our studies and the re- 
sults lead basically to the same general conclusions, albeit 
some important differences should be noted. In our analy- 
sis of the index distribution within a single burst (Fig. 0), 
the rji — rjj differences show a narr ower distribution than 
the corresponding 7$ — 7^ . In § |4.l| we concluded that the 
new approach gives better correlations. This implies that, 
in general, the estimated uncertainties of r\ are smaller. 
Even so, when testingthc null hypothesis of index invari- 
ance using equation (U3) , Ho is marginally accepted in the 
7 case (probability p = 7%) , while it is fully consistent for 
77 (p = 64%). We have also shown that the general distri- 
butions of 7/ and 7 (Fig. H) are very similar. Their means 
and dispersions are equal to within the uncertainties, thus 
covering about the same range of values. All these facts 
strongly suggest that the underlying distribution within a 
single burst is indeed very narrow, that the observed larger 
spread can be attributed to the uncertainties of the mea- 
surements, and that a better estimator of the HIC will 
consequently reduce this spread. Obviously, this is also 
valid for the determination of the general distribution, but 
in this case the observed reduction is comparable to the 
standard deviation uncert ainti es. 

We also discussed in 8 4.4 that within a burst, later 



comparative simulations of the methods, assuminga con- 
stant spectral shape (i.e., a and /3 fixed in eq. |g|). A 
larger sample is needed to verify the trend, but this sug- 
gests that indeed the new method is less affected by the 
detector energy window. 

The tp-m.eih.od thus provides a better way of studying 
detailed features in the observed HICs. Finally, the better 
invariance within a burst shown by the 77 indices may al- 
low us to recognize hidden pulses in the light curve of the 
burst, as will be discussed next. 

5.2. Track Jumps and Overlapping Pulses 

Figure [}p illustrates the observed evolution of (p(E p ^) 
along a track during the decay phase of a single pulse. 
When many pulses are present in the light curve of a burst, 
and even though there is no apparent general trend dur- 
ing the rise phases, we have found that they follow almost 
parallel tracks during the decay phases. To observe this be- 
havior clearly, it is necessary that the pulses do not overlap 
significantly. Also, they must be bright and long enough 
to establish t he dir ection of the track, see section 3.2. 

In section [f.3.1 , we presented the case GRB 950104 
(trigger 3345) with marginally separated pulses, which ex- 
hibits a track jump in its HIC. We consider this as a limit- 
ing case, where it can still be claimed, from the analysis of 
the time history alone, that a secondary pulse is present 
(see Fig. ga). This fact, together with the slope agreement, 
within the uncertainties, between the first pulse track and 
the tracks in the second pulse, makes it an interesting ex- 
ample among the track jump cases. 

Analogous HIC behavi ors we re found in a further three 
cases, shown in section 4.3.2, but there no visible sec- 



pulses seem to have larger 7 indices, although the differ- 
ence is not highly significant. This trend is not observed 
in the 77 values. This might be a systematic error over- 
come by the c/j-method. Later pulses are usually softer, 
and therefore closer to the lower energy limit of observa- 
tion. See, for example, the burst trigger 2083 in Figure [I[ 
The greater slope of the second pulse PLHIC, quantified 
in panel (c), is of the same order as the one obtained in 



ondary pulse was seen. For instance, GRB 970925 (trig- 
ger 6397; see Fig. [3p) shows an apparent change in the 
overall decay, which could be due to a hidden pulse. How- 
ever, there is an implicit assumption of the pulse shape 
being an exponential. Nevertheless, in the present study 
we do not want to make such assumptions concerning the 
shape of the pulses. 

A simple interpretation of the track jump is to assume 
that it is produced by an overlapping secondary pulse. 
From this point of view, this feature is present also in cases 
like GRB 970420 (trigger 6198) shown in Figure |. The 
only difference is that there we can resolve the individual 
pulses in the light curve. Due to the strong temporal evo- 
lution of -Epk observed in all presumed single pulses, two 
overlapping pulses will show, if all pulses within a burst 
have similar spectral shapes and behaviors, a continuous 
change in the way the individual spectral shapes overlap 
each other. These fundamental shapes may be fairly con- 
stant, but the combined evolution in the EF^-E plane 
will give the impression of an overall change. In Figure ph, 
we illustrate such a model of two spectral components. 
The components follow identical tracks during their decay 
phases. As the peak energy of the main pulse declines, 
a harder component appears at the rise phase of the sec- 
ondary pulse. At some point this component starts to 
dominate the total spectrum and the fitting routine finds 
a better x 2 value, shifting the Band et al. function param- 
eters to peak at a higher E p k, and the jump occurs. Note 
that due to the broad shape of the spectrum, the changes 
in the values of the asymptotic slopes of the Band et al. 
function are an artifact of the fitting. 
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The uncertainties of the parameter £" p k are estimated 
by evaluating the confidence region around the \ 2 minima 
of the fits. When the spectral shape has a local maxi- 
mum close to the absolute maximum, an overestimation 
of the -Epk uncertainties occurs. These uncertainties are of 
the same order as the energy interval that separates the 
tracks. But the very good alignment of the data points 
along the tracks indicates that the shape features are real 
and not a random occurrence. 

In studies of the spectral/temporal evolution, other 
workers (e.g., Kargatis et al. 1995; Liang & Kargatis 1996) 
have fixed the parameters a and p to reduce the estimated 
uncertainties of -Epk- For instance, the parameter values 
of the time-integrated spectrum can be used. Such a so- 
lution provides more "stable" -E p k values, but a case like 
GRB 950104 (trigger 3345) would be missed. It smoothes 
out the track jump in such a way that a single, steeper 
track is obtained for the whole decay. Also a coarser time 
binning (see Ryde & Svensson 2000c) makes the jump "in- 
visible" . 

Based on this discussion we propose a pulse identifica- 
tion method for cases of extensive overlap, under the fol- 
lowing assumptions. First, the GRB pulses are assumed to 
actually follow a PLHIC during their decays. We may not 
observe this because the pulse is too weak or short, because 
the observation energy window is inadequate, or because 
there are unresolved secondary pulses present. Second, we 
have to assume that the correlation index is approximately 
constant within a burst. We have shown four cases where 
the pulse separation is evident in the hardness-intensity 
plane but not as clearly in the light curve. This pulse sepa- 
ration's not so easily accomplished by traditional methods 
(see §3.2), but can be used as an auxiliary tool. However, 
to identify track jumps with some confidence, the distance 
between tracks should be larger than the observed disper- 
sion along individual tracks. Also a minimum number of 
points on each track is needed to produce a reasonable fit. 

We found few cases in our BATSE sample that meet 
these requirements. However, future instrumentation will 
allow higher time-spectral resolution and these criteria 
may then have a wider application. 

6. CONSTRAINTS ON GRB MODELS 

The conditions that determine a particular PLHIC index 
may depend, in principle, on many physical factors of the 
emitting system, e.g., the density, composition, radiation 
processes, geometry, etc. This fact is reflected by the broad 
distribution of the indices r\ and 7 that we found when 
studying pulses from dif ferent bursts (see Fig . |3|) . The in 
tcrnal models (see , e.g 



Kobayashi et al. 1997, Daigne & 



Mochkovitch 199?) rely on the assumption of a random 
distribution of some shell parameters, such as the den- 
sity, the kinetic energy or the momentum, to explain the 
chaotic time histories observed but also to obtain high ef- 
ficiencies. The sequence of collisions and merging between 
different shells follows a stochastic process that will hardly 
consistently reproduce identical physical conditions in each 
emitting shock front. On the contrary, most probably it 
will statistically reflect almost the same diversity of condi- 
tions found when we observe many bursts. Based on these 
models, one would expect a similar parameter dispersion 
when pulses are taken from the same or different bursts. 
On the other hand, it is conceivable that in the case bursts 



were very anisotropic, differences due to the viewing angle 
would be enhanced by the relativistic beaming and this 
could become a dominating factor. In spite of that, the 
internal models will still need to produce a narrow PL- 
HIC index dispersion within a burst. The approximate 
constancy of the PLHIC index is more easily explained in 
terms of a single system where an active region is regen- 
erated or different parts become active at different times, 
as it is assumed in many variations of t he external shock 
scenario (see, e.g., Fenimore et al. 1999| ). 



7. CONCLUSIONS 

We have introduced a new method to study the 
hardness- intensity correlation in GRBs. Instead of mea- 
suring the energy flux, we use the value of ip = EF^, at 
E = -Epk, which differs by an approximately constant fac- 
tor. This method has many advantages. First of all, the 
correlations become stronger. Second, the tp-method is 
less dependent on the spectral window of the observations. 
Third, it disregards the effects of soft components and ad- 
ditional weak pulses with spectra which overlap the one 
being studied. In cases where the spectrum consists of 
two approximately equally strong components (from dif- 
ferent pulses) the ip-method captures the evolution of one 
of the components and disregards the second. These cases 
produce broad spectra which are not very well modeled by 
the Band et al. function. Fourth, the resulting flux mea- 
sure, f, is less sensitive to the details of the deconvolution, 
e.g., the uncertainties in the spectral shape. 

We constructed a sample of prominent GRB pulses con- 
sisting of 82 cases. Of these, we found that at least 75% 
exhibited moderately good (or better) PLHICs. Among 
the poorly correlated cases, we could recognize some tech- 
nical problems that were affecting our analysis. Thus, the 
applicability of the PLHIC model may be even larger. 

A specific feature was found in the hardness-intensity 
evolution of some seemingly single pulses. In the HIC di- 
agram of these cases, there is a sharp transition from one 
power law track to a second, both with the same PLHIC 
index. This feature is denoted as a track jump and in half 
of the cases, it coincides with a weak feature in the light 
curve. This could be explained as being produced by over- 
lapping pulses. The track jumps in the HIC diagram reveal 
the transition to a new pulse and it can be used as an aid 
for pulse identification. 

We established that the distribution of the PLHIC index 
is narrower for pulses within a burst than the distribution 
for pulses from different bursts. The latter can be approx- 
imated by a normal distribution. The dispersion within a 
burst could be entirely attributed to the stochastic uncer- 
tainties and the index is thus practically invariant. 

These results demand a physical model to be able to 
reproduce multiple pulses from an individual burst with 
similar characteristics. At the same time, however, pulses 
from an ensemble of different bursts should exhibit a di- 
versity that is larger. This is particularly relevant when 
comparing the external versus the internal models. The 
latter models require a large diversity of properties of the 
colliding shells giving rise to the pulses, both within a burst 
and among bursts. 

Changes in the observational conditions among differ- 
ent bursts, such as viewing angle, could possibly account 
for the disparity between the index distributions. These 
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matters will be considered in future research. 
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Fig. 1. — Examples of pulses from our sample. The upper two rows show fairly smooth, typical pulses selected. However, we also include 
some cases that show some apparent sub-structure, as the ones in the bottom row. See Table hi for further information. 




Fig. 2 r — Relation between the correlation indices 7 and 77 for 47 pulses with a good PLHIC, i.e., A77/77 < 0.15 (marked as subset A 
in Table El). The straight line shows a weighted linear fit 7 = (0.96 ± 0.09)»; + (0.08 ± 0.16) to the data. This shows that there is a good 
correspondence between the two methods, i.e., using the intensity measures F and ip, respectively. 
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Fig. 3. — Distributions of the hardness-intensity correlation indices, 7 and r;, for a set of 39 pulses (marked in Table U as subset B). This 
set is defined by all the sample cases with Arj/r; < 0.15, but restricted to only one pulse per burst (see text for more (retails). Both index 



distributions can be approximated by normal distributions with 7 
respectively. 



1.9 ± 0.1 and <x 7 = 0.72 ± 0.08, and fj = 2.0 ± 0.1 and cr„ = 0.68 ± 0.08 
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Fig. 4. — a) Count rate time history showing the time binning used to study the spectral evolution of GRB 921207 (trigger 2083). A few bin 
numbers are shown. The last few bins are not included in this study because the _E p i5 estimation is very uncertain as the peak of the spectrum 
approaches the low energy limit of the detector. This is often seen when i? pk < 40 keV. b) Time evolution of the ratio F/ip. The integration 
of the flux was made over the energy range 24—2000 keV. The F'-measure will suffer from window problems at high and low energies, c) The 
PLHIC measured with the </j-method. The evolution during the decay phases is shown with dark dots. Error bars are excluded for clarity, and 
only a typical average value is shown separately. The numbers refer to the time- bins in panel (a). We find, for each pulse decay, the power 
law indices rji = 1.78 ± 0.04 and r\2 = 1.77 ± 0.06. d) Same as (c), now using the F-method. The corresponding indices are 71 = 1.74 ± 0.04 
and 72 = 2.03 ± 0.07. Note that, as in general, the similarity between the PLHIC indices of pulses within a burst is more evident using the 
</3-method. 
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Fig. 5. — a) Count time history of GRB 970420 (trigger 6198) showing the binning used to study the spectral evolution. The second pulse 
(at bin number 8) actually consists of two close-lying pulses, which are seen in the 64-ms data, b) The peaks of the instantaneous spectra at 
successive times are represented by dots. Error bars are excluded for better visibility, and only a typical average value is shown separately. 
The numbers refer to the time- bins in panel (a). During the decay phases (dark dots), the HIC follows very good power laws with indices 
r\\ = 2.25 ± 0.20, r]2 = 2.40 ± 0.15, and r?3 = 2.20 ± 0.15, equal to within the estimated uncertainties. 
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GRB 950104 (trigger 3345) 
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Fig. 6. — Spectral evolution of GRB 950104 (trigger 3345) and GRB 970925 (trigger 6397). a), d) Photon time history showing the binning 
used to study the spectral evolution (S/N = 40). The dashed line indicates the background level, b) The same time interval using 64-ms time 
resolution data. A secondary peak is present at t ~ 6.5s. See the text for details, c) The peaks of the instantaneous spectra at successive 
times are represented by dots. The numbers refer to the time-bins. The first pulse decay follows a power law with index rji = 0.77± 0.06. The 
second pulse decay initially follows the same law, within the measurement uncertainties, with r)2 = 0.78 ± 0.025. A track jump occurs at the 
time of the rise of the secondary peak in the light curve and then the spectral peak evolution continues on a parallel track rjs = 0.81 ± 0.05 
while the light curve decays again. Error bars are excluded for better visibility, and only a typical average value is shown separately. The 
fitting procedure largely overestimates the uncertainties in these cases, e) The burst trigger 6397 at 64-ms time resolution. /) Evolution of the 
peak energy during the pulse decay phase. As in the previous case (c), the whole decay presents a poor hardness-intensity correlation, but it 
can be divided into parallel tracks with very good HICs and equal indices to within the uncertainties (rji = 0.80 ± 0.15 and 772 = 0.90 ± 0.07). 
Note the presence of a small "bump" at t ~ 7.5 s, i.e., at bins number 6 and 7. 



1G 



The HIC in GRB pulses 



GRB 910927 (trigger 829) 




GRB 960912 (trigger 5601) 



time [s] 



1000 


b> if \ 


64-ms 


800 




; 


600 




: 


400 




; 


200 




yyyxhmiih 



2.5 5 7.5 10 12.5 15 17.5 20 

time [s] 




17.5 20 



700 


c) 


/ / 2 
9/ // 1 


500 




L**f' i 


(N 






E 






o 




12/ / 


V 300 






w 






> 






<D 






^-200 






LU 




/ 14 


u_ 






lu 150 






100 


/ 7 16 


<AEpeak> = * 





f) 






/*9 


500 






12/' 




■\T^ 






/• /^ 




§ 200 






-/ ' /^ 




'co 






^/18 




> 

J[100 










LU 
LL 










m 50 




18 






20 






<AEpeak>= . — 





70 100 

Epeak [keV] 



150 



200 



50 100 200 

Epeak [keV] 



500 



1000 



Fig. 7.— Spectral evolution of GRB 910927 (trigger 829) and GRB 960912 (trigger 5601). a), d) Photon time history showing the binning 
used to study the spectral evolution with signal-to-noise levels of 40 and 30, respectively. The dashed line indicates the background level. 
6), e) Time history at 64-ms time resolution, c) Evolution of the peak energy during the pulse decay phase. As in the previous examples, 
it can be divided into two parallel tracks with indices r)\ = 2.95 ± 0.3 and r)2 = 2.91 ± 0.07. Note that the first pulse, that has its peak at 
bin 2, seems to follow the track defined by the last three data points numbered 14-16. It could therefore be that the pulse peaking at bin 9 
has a fast decline and is amid the decay of a dominating, long pulse. /) The same for burst trigger 5601, with indices r]i = 0.675 ± 0.01 and 
r\2 = 0.675 ± 0.04, equal to within the uncertainties. In both cases the high-resolution light curve shows no clear sign of a "hidden" pulse. 
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error bars of the r\ differences are slightly smaller than those of 7. Furthermore, most error bars cover the origin, implying equal values to 
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Table 1 
Sample of 82 pulses in 66 GRBs 



Burst 


trigger 


LAD 


WxH 


"bins 


GRB910627 


451 


4 


5.6 


5 


GRB910897 


647 





3.2/14.0 


8/8 


GRB910814 


676 


2 


51.5 


5 


GRB910927 


829 


4 


6.3 


8 T.T 


GRB911016 


907 


1 


1.4 


14 


GRB911031 


973 


3 


2.7/24.2 


19/12 


GRB911104 


999 


2 


4.0 


5 


GRB911118 


1085 


4 


8.7 


14 


GRB911126 


1121 


4 


21.9/24.7 


8/6 


GRB911202 


1141 


7 


4.0/9.7 


8/8 


GRB920221 


1425 


7 


5.9 


5 


GRB920502 


1578 


5 


8.0 


9 


GRB920525 


1625 


4 


5.0/12.6 


12/5 


GRB920623 


1663 


4 


17.1 


10 


GRB920718 


1709 


7 


2.1 


6 


GRB920830 


1883 





0.8 


7 


GRB920902 


1886 


5 


8.3 


10 


GRB921003 


1974 


2 


1.3 


5 


GRB921123 


2067 


1 


21.2 


8 


GRB921207 


2083 





1.1/8.6 


22/9 


GRB930201 


2156 


1 


14.9 


12 


GRB930425 


2316 


1 


14.4 


12 


GRB930612 


2387 


6 


5.9 


19 


GRB931221 


2700 


3 


53.4 


8 


GRB940329 


2897 


4 


8.2/24.0 


5/5 


GRB940410 


2919 





0.5 


13 


GRB940429 


2953 


3 


5.7/13.0 


9/7 


GRB940623 


3042 


1 


7.8 


14 


GRB940708 


3067 





2.4 


19 


GRB941026 


3257 





4.5 


13 


GRB941121 


3290 





39.7 


10 


GRB950104 


3345 


1 


5.8 


-. T.T 


GRB950325 


3480 


3 


0.2 


10 


GRB950403 


3491 


3 


7.7 


16 


GRB950403 


3492 


5 


5.3 


15 


GRB950624 


3648 


3 


2.7/22.7/40.9 


6/5/10 


GRB950818 


3765 


1 


66.1 


12 


GRB950909 


3788 


3 


27.2 


14 


GRB951016 


3870 


5 


0.5 


8 


GRB951102 


3891 


2 


33.3 


7 


GRB951213 


3954 


2 


0.8 


16 


GRB960113 


4350 


1 


13.8 


6 


GRB960124 


4556 


5 


1.8/3.6 


6/5 


GRB960530 


5478 


2 


1.9 


7 


GRB960531 


5479 





52.5 


5 


GRB960708 


5534 


5 


1.7 


5 


GRB960804 


5563 


4 


1.4 


6 


GRB960807 


5567 





11.8 


10 


GRB960912 


5601 





1.9 


10 TJ 


GRB961001 


5621 


2 


3.9/7.2 


6/5 


GRB961009 


5628 


1 


10.2 


5 


GRB961102 


5654 


5 


19.0 


14 


GRB961126 


5697 


6 


0.6 


6 


GRB970111 


5773 





8.1/17.3/19.4 


8/5/6 


GRB970223 


6100 


6 


8.2 


16 


GRB970420 


6198 


4 


4.2/6.1/9.9 


5/8/5 


GRB970815 


6335 


7 


0.8 


7 


GRB970925 


6397 


7 


2.8 


12 TJ 


GRB971127 


6504 


2 


3.2 


5 


GRB980125 


6581 





47.6 


7 


GRB980301 


6621 


1 


32.7 


(i 


GRB980306 


6629 


1 


215 


(> 


GRB980306 


6630 


3 


2.1 


9 


GRB980821 


7012 





2.9 


6 


GRB990102 


7293 


6 


3.3 


10 


GRB990123 


7343 





37.7 


8 



Case with a track jump feature in its hardness-intensity tem- 
poral evolution 



Borgonovo & Ryde 



19 



Table 2 
The HIC power law indices for 78 pulses in 62 GRBs 



Pulse 



7 



Hi 



II: 



Pn (-R ) Arj/rj Comments 



451b 


1.36 ±0.24 


1.38 ±0.33 


0.9166 


0.8509 


647a 


2.12 ± 0.22 


1.80 ±0.25 


0.9382 


0.8868 


647b 


2.54 ±0.18 


1.40 ±0.30 


0.9692 


0.7888 


676 


3.30 ±0.40 


2.80 ±0.25 


0.9561 


0.9788 


907 


2.34 ±0.15 


2.54 ±0.20 


0.9512 


0.9422 


973a 


1.45 ± 0.21 


1.28 ±0.21 


0.7371 


0.6821 


973b 


0.84 ± 0.10 


0.71 ±0.11 


0.8590 


0.8041 


999 


1.90 ±0.10 


1.60 ±0.10 


0.9905 


0.9825 


1085 


1.79 ±0.025 


2.05 ±0.04 


0.9978 


0.9954 


1121a 


1.30 ±0.15 


1.00 ±0.10 


0.9318 


0.9172 


1121b 


1.17 ±0.20 


1.10 ±0.14 


0.9150 


0.9385 


1141a 


0.84 ±0.10 


0.51 ±0.10 


0.9198 


0.8117 


1141b 


0.86 ±0.15 


0.76 ±0.15 


0.8781 


0.8400 


1425 


2.70 ±0.50 


2.55 ±0.50 


0.8922 


0.9060 


1578 


2.80 ±0.40 


2.95 ±0.55 


0.8696 


0.8038 


1625a 


2.37 ±0.30 


2.12 ±0.31 


0.8584 


0.8188 


1625b 


2.23 ±0.58 


1.95 ±0.55 


0.8315 


0.8057 


1663 


1.45 ±0.15 


1.00 ±0.15 


0.9139 


0.8342 


1709 


2.30 ±0.15 


2.32 ±0.08 


0.9865 


0.9946 


1883 


2.39 ±0.29 


2.51 ±0.40 


0.9303 


0.8783 


1886 


3.80 ±0.45 


3.60 ±0.45 


0.8965 


0.8786 


1974a 


2.20 ±0.30 


2.10 ±0.30 


0.9537 


0.9495 


2067 


2.60 ±0.15 


2.60 ±0.10 


0.9852 


0.9902 


2083a 


1.78 ±0.04 


1.74 ±0.04 


0.9889 


0.9890 


2083b 


1.77 ±0.06 


2.03 ±0.07 


0.9915 


0.9926 


2156 


1.55 ±0.15 


1.00 ±0.13 


0.9117 


0.8460 


2316 


0.64 ±0.23 


0.60 ±0.30 


0.4395 


0.2870 


2387 


1.09 ±0.17 


0.83 ±0.33 


0.7035 


0.2695 


2700 


1.70 ±0.45 


1.70 ±0.45 


0.7067 


0.6867 


2897b 


1.24 ± 0.45 


1.47 ±0.40 


0.7344 


0.7886 


2897c 


1.42 ±0.40 


2.20 ±0.50 


0.8177 


0.8179 


2919 


1.39 ±0.38 


2.20 ±0.50 


0.5478 


0.3954 


2953a 


1.61 ±0.30 


1.28 ±0.27 


0.8094 


0.7577 


2953b 


1.46 ±0.45 


1.24 ±0.41 


0.6785 


0.5777 


3042 


1.38 ±0.88 


1.07 ±0.90 


0.3809 


0.7875 


3067 


1.35 ±0.10 


1.47 ±0.13 


0.8963 


0.8820 


3257 


2.80 ±0.20 


2.60 ±0.20 


0.9573 


0.9268 


3290 


1.12 ±0.08 


1.41 ±0.10 


0.9625 


0.9560 


3480 


2.07 ±0.15 


1.90 ±0.13 


0.9686 


0.9653 


3491 


1.75 ± 0.09 


1.59 ±0.10 


0.9600 


0.9477 


3492 


1.65 ±0.10 


1.50 ±0.10 


0.9378 


0.9272 


3648a 


1.03 ±0.30 


1.87 ±0.50 


0.7705 


0.7647 


3648b 


0.64 ±0.03 


1.24 ±0.25 


0.9930 


0.9027 


3648c 


1.42 ±0.08 


1.38 ±0.08 


0.9743 


0.9749 


3765 


2.40 ± 0.20 


2.45 ±0.20 


0.9438 


0.9217 


3788 


1.67 ±0.26 


1.35 ±0.35 


0.7691 


0.5571 


3870 


0.78 ±0.15 


0.76 ±0.20 


0.7927 


0.7147 


3891 


1.80 ±0.18 


1.58 ±0.25 


0.9498 


0.8922 


3954 


1.25 ± 0.20 


0.97 ±0.21 


0.7759 


0.5907 


4350 


0.82 ±2.8 


1.0 ±4.3 


0.0207 


0.0137 


4556a 


1.50 ±0.30 


1.20 ±0.30 


0.8756 


0.7937 


4556b 


1.47 ±0.17 


1.36 ±0.23 


0.9621 


0.9180 


5478 


1.87 ±0.16 


1.05 ±0.20 


0.9664 


0.8835 


5479 


1.06 ±0.26 


1.03 ±0.30 


0.8423 


0.7785 


5534 


0.55 ±0.14 


0.96 ±0.17 


0.8415 


0.9138 


5563 


2.35 ±0.30 


2.40 ±0.40 


0.9226 


0.9042 


5567 


3.96 ±0.53 


3.68 ±0.50 


0.8750 


0.8761 


5621a 


1.26 ±0.26 


0.97 ±0.32 


0.8521 


0.6893 


5621b 


1.32 ±0.16 


1.35 ±0.12 


0.9570 


0.9773 


5628 


1.05 ±0.25 


0.75 ±0.20 


0.8503 


0.7850 


5654 


2.00 ±0.50 


2.00 ±0.50 


0.5549 


0.5728 


5697 


3.75 ± 0.70 


4.50 ±0.70 


0.8804 


0.9057 


5773a 


1.56 ±0.15 


1.35 ±0.14 


0.9573 


0.9397 


5773b 


1.82 ±0.10 


1.32 ±0.50 


0.9868 


0.7122 


5773c 


1.58 ±0.15 


1.60 ±0.20 


0.9614 


0.9261 


6100 


2.03 ±0.18 


1.79 ±0.20 


0.8992 


0.8549 


6198a 


2.25 ± 0.20 


2.01 ±0.18 


0.9793 


0.9769 


6198b 


2.38 ±0.15 


2.16 ±0.14 


0.9784 


0.9766 


6198c 


2.18 ± 0.12 


2.29 ±0.14 


0.9908 


0.9892 


6335a 


2.35 ±0.20 


2.70 ±0.30 


0.9727 


0.9431 


6504 


1.96 ±0.36 


2.38 ±0.70 


0.9068 


0.8117 


6581 


1.62 ±0.20 


1.40 ±0.20 


0.9374 


0.9327 


6621 


2.0 ± 1.0 


2.30 ± 1.2 


0.4892 


0.4738 


6629 


1.20 ±0.9 


0.95 ±0.85 


0.3274 


0.2319 


6630 


1.49 ±0.04 


1.58 ±0.05 


0.9944 


0.9935 


7012 


1.91 ± 0.44 


1.60 ±0.40 


0.8231 


0.7448 


7293 


2.04 ± 0.12 


2.05 ±0.22 


0.9731 


0.9145 



1.0 x IO" 2 
7.6 x IO" 5 

9.2 x IO" 6 

4.0 x 10" 3 

3.1 x 10" 9 

2.6 x 10" 6 

1.5 x 10~ 5 
3.9 x 10" 4 

< io- 10 

1.0 x IO" 4 

2.8 x IO" 3 

1.7 x IO" 4 

5.9 x IO" 4 

1.6 x IO" 2 
2.5 x IO" 4 

1.5 x IO" 5 

3.1 x IO" 2 

1.6 x IO" 5 
6.9 x IO" 5 
4.5 x IO" 4 

3.3 x IO" 5 

4.3 x IO" 3 
1.0 x IO" 6 

< io- 10 

1.7 x IO" 8 

1.4 x IO" 6 

1.8 x IO" 2 

7.2 x IO" 6 

8.9 x IO" 3 

6.4 x IO" 2 

3.5 x IO" 2 

3.8 x IO" 3 

9.5 x IO" 4 

2.3 x IO" 2 

1.9 x IO" 1 

8.6 x IO" 10 

7.0 x IO" 9 
5.5 x IO" 7 

2.7 x IO" 7 

< IO" 10 
3.2 x IO" 9 

2.1 x IO" 2 

2.5 x IO" 4 

1.2 x IO" 7 

1.4 x IO" 7 

3.8 x IO" 5 
3.0 x IO" 3 

2.0 x IO" 4 

6.6 x IO" 6 

7.9 x IO" 1 

6.1 x IO" 3 

3.2 x IO" 3 

7.1 x IO" 5 
2.8 x IO" 2 
2.8 x IO" 2 

2.3 x IO" 3 
7.0 x IO" 5 
8.6 x IO" 3 
3.8 x IO" 3 
2.6 x IO" 2 

2.2 x IO" 3 

5.6 x IO" 3 

2.5 x IO" 5 
6.5 x IO" 4 

5.7 x IO" 4 

2.3 x IO" 8 

1.3 x IO" 3 
3.2 x IO" 6 

3.8 x IO" 4 
4.2 x IO" 5 
1.2 x IO" 2 

3.4 x IO" 4 
1.2 x IO" 1 
2.4 x IO" 1 
3.8 x IO" 9 
1.2 x IO" 2 
1.4 x IO" 7 



0.18 




0.10 


A 


0.071 


A,B 


0.12 


A,B 


0.064 


A,B 


0.14 


A,B 


0.12 


A 


0.053 


A,B 


0.014 


A,B 


0.12 


A,B 


0.17 




0.12 


A,B 


0.17 




0.19 




0.14 


A,B 


0.13 


A,B 


0.26 




0.10 


A,B 


0.065 


A,B 


0.12 


A,B 


0.12 


A,B 


0.14 


A,B 


0.058 


A,B 


0.022 


A,B 


0.034 


A 


0.097 


A,B 


0.35 




0.16 




0.26 




0.36 




0.28 




0.27 




0.19 




0.31 




0.64 




0.074 


A,B 


0.071 


A,B 


0.071 


A,B 


0.072 


A,B 


0.051 


A,B 


0.061 


A,B 


0.29 




0.047 


A 


0.056 


A,B 


0.083 


A,B 


0.16 




0.19 




0.10 


A,B 


0.16 


(1) 


3.4 




0.20 




0.12 


A,B 


0.086 


A,B 


0.25 




0.25 




0.13 


A,B 


0.13 


A,B 


0.21 




0.12 


A,B 


0.24 




0.25 




0.19 




0.096 


A,B 


0.055 


A 


0.095 


A 


0.089 


A,B 


0.089 


A 


0.063 


A,B 


0.055 


A 


0.085 


A,B 


0.18 




0.12 


A,B 


0.50 




0.75 




0.027 


A,B 


0.23 




0.059 


A,B 



2(1 



The HIC in GRB pulses 



Table 2 — Continued 



Pulse 



r(R 2 v ) 



Ari/rj 



Comments 



7:] 13 



1.50 ±0.50 



1.10 ±0.50 



0.5905 



0.4201 



0.33 



Note. — Pulses are denoted by the burst trigger number and, when it is necessary, labeled alphabetically following 
temporal order. The labeling is consistent with the additional pulses presented in Table |3J. In the comments column 
we specify the pulses belonging to the subset A consisting of 47 pulses from 39 bursts, and defined by the condition 
Ar;/?? < 0.15. The subset B is also indicated, and it is constructed as A, but taking only one pulse per burst (choosing 
the onii with lowest Pjv(-R^)), therefore it comprises 39 pulses from 39 bursts. Note that the four track jump cases from 
Table III arc excluded here. 

Rejected as a is not constrained 



Table 3 

Extended Sample of Multi-Pulse Bursts 



Pulse 


WW 


"bins 


V 


7 


R 2 


Pn{R%) 


Ari/r/ 


451a 


0.0 


4 


0.98 ±0.04 


0.70 ±0.12 


0.9963 


1.9 x 10~ 3 


0.041 


451b 


5.6 


5 


1.36 ±0.24 


1.38 ±0.33 


0.9166 


1.0 x 10~ 2 


0.18 


1974a 


1.3 


5 


2.20 ±0.30 


2.10 ±0.30 


0.9537 


4.3 x 10 -3 


0.14 


1974b 


6.5 


4 


2.15 ±0.55 


3.05 ±0.70 


0.8796 


6.2 x 10~ 2 


0.26 


2897a 


5.9 


4 


1.13 ±0.40 


1.30 ±0.40 


0.7062 


1.6 x 10 _1 


0.35 


2897b 


8.2 


5 


1.24 ±0.45 


1.47 ±0.40 


0.7344 


6.4 x 10~ 2 


0.36 


2897c 


24.0 


5 


1.42 ±0.40 


2.20 ±0.50 


0.8177 


3.5 x 10~ 2 


0.28 


3345a 


3.5 


4 


0.77 ±0.05 


0.60 ±0.06 


0.9909 


4.6 x 10 -3 


0.065 


3345b 


5.8 


3 


0.79 ±0.03 


0.61 ± 0.065 


0.9985 


2.5 x 10~ 2 


0.038 


3345c 


6.5 


4 


0.76 ±0.05 


0.62 ±0.17 


0.9892 


5.4 x 10~ 3 


0.066 



NOTE. — The multi-pulse sample comprises 35 pulses from 15 bursts, and it is an extension of the main sample that 
includes also 4 time-bin pulses (not listed in Table 111). All the bursts studied in this set have at least one pulse having more 
than 4 time-bins that belong to the main sample. Tlcrc we present only the bursts with additional pulses for conciseness. 
The whole set is, listed by burst trigger number (number of pulses): 451(2) , 647(2), 973(2), 1121(2), 1141(2), 1625(e), 
1974(2), 2083(2), 2S97(3), 3345(3), 3648(3), 4556(2), 5621(2), 5773(3), 6198(3) (cases in italic appeared in Table g|) . 
Pulse 3345 in Table 111 is listed as a track jump case, and it has been separated into two short pulses (b and c) with 3 and 
4 timc-bins, respectively (see the discussion in the text). 



Table 4 

Pulses with Track Jumps 



Pulse 


V 


7 


*; 


Pn{R%) 


Ari/r, 


829a 


2.87 ±0.30 


3.20 ±0.12 


0.9816 


9.2 x 10~~ 3 


0.10 


829b 


2.91 ±0.06 


3.50 ±0.25 


0.9995 


2.5 x 10~ 4 


0.021 


5601a 


0.67 ±0.013 


0.48 ±0.07 


0.9992 


4.0 x 10~ 4 


0.019 


5601b 


0.68 ±0.04 


0.74 ±0.03 


0.9966 


3.7 x 10~ 2 


0.059 


6397a 


0.77 ±0.13 


0.24 ±0.25 


0.9265 


8.7 x 10~ 3 


0.17 


6397b 


0.86 ±0.05 


0.49 ±0.06 


0.9889 


5.0 x 10~ 4 


0.058 



Table 5 
Hypothesis Tests - P- values 



Test 


1 


7 


Xl (Eq. bj) 
T a (Eq. M) 
Xl (Eq. M) 
T b (Eq. Q) 


10- 7 
< 10- 5 

7 x 10~ 14 ||0.38 * 
< 10~ 5 ||0.64 * 


5 X 10~ 6 
4 x 10~ 5 

0.03 

0.07 



* Results with and without the outlier data 
point discussed in the text 



